Dynamics of a Brownian circle swimmer 
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Self-propelled particles move along circles rather than along a straight line when their driving 
force does not coincide with their propagation direction. Examples include confined bacteria and 
spermatozoa, catalytically driven nanorods, active, anisotropic colloidal particles and vibrated gran- 
ulates. Using a non-Hamiltonian rate theory and computer simulations, we study the motion of a 
Brownian "circle swimmer" in a confining channel. A sliding mode close to the wall leads to a huge 
acceleration as compared to the bulk motion, which can further be enhanced by an optimal effective 
torque-to-force ratio. 

PACS numbers: 05.40.Jc, 82.70.Dd 



Active particles, which are self-propelled by their own 
motor, exhibit a wealth of novel and fascinating nonequi- 
librium effects such as giant density fluctuations 
swarming and swirling 3]. Examples are found 
in quite different areas of physics and include micro- 
organisms propelled by flagella in a fluid 0, H, H, 5 H] > 
man-made colloidal swimmers J9fl, catalytically driven 
nanorods or Janus particles fiol . Illj | , vibrated granulates 
of polar rods [!, [l| , and pedestrians [l3[ . Typically it is 
assumed that the swimmers move along their symmetry 
axis such that the force and the particle orientation are 
in line. This leads to a motion along a straight line just 
perturbed by random (e.g., Brownian) fluctuations. 

Here we study the case in which the internal force pro- 
pelling a colloidal particle does not coincide with the 
particle orientation. In the absence of Brownian fluc- 
tuations, this will lead to an overdamped motion along a 
closed circle, therefore we refer to this particle as a "cir- 
cle swimmer." Even a slight misalignment of the drive 
direction will result in circle swimming, which is thus the 
generic case of self-propulsion. Circle swimmers with a 
pronounced curved trajectory are realized in nature and 
can be artificially prepared: In fact, it has been shown 
that certain bacteria [g, 0, H, EH and spermatozoa 0, @ , 
when confined to two dimensions, swim in circles. More- 
over, catalytically driven nanorods [1(3, EH and colloidal 
particles [9( can be prepared with a tilted motor, and 
a vibrated polar rod 0] on a planar substrate with an 
additional left-right asymmetry will move along circles. 
Last but not least, the trajectories of completely blinded 
and ear-plugged pedestrians have a significant circular 
form [l5|]. Despite their practical importance, the Brow- 
nian dynamics of a circle swimmer has not yet been ad- 
dressed by theory and simulation either in the bulk or 
under confinement flil ]. 

In this paper, we propose a simple model for Brown- 
ian motion of a circle swimmer in two spatial dimensions 
arising from the combined actions of an internal self- 
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propelling force and a torque. We solve the Langevin 
equation of a two-dimensional circle swimmer analyti- 
cally in the bulk providing a suitable reference model. 
The averaged position falls on a spira mirabilis, and a 
crossover from an oscillatory ballistic to a diffusive be- 
havior is found in the mean-squared displacement. We 
then identify the modes of propagation of a circle swim- 
mer in confining channels with repulsive walls using com- 
puter simulations and a non-Hamiltonian rate theory. In 
symmetric channels, the long-time self-diffusion coeffi- 
cient Dl is significantly enhanced mediated by an effi- 
cient sliding mode of a tilted rod close to a wall. Fur- 
thermore, Dt, is nonmonotonic in the torque. Finally, in 
asymmetric channels which are lacking a left-right sym- 
metry (e.g., due to gravity [l7j|). the sliding mode of the 
circle swimmer yields a ballistic motion along the wall. 

Neglecting hydrodynamic interactions, the over- 
damped motion of the Brownian circle swimmer in two 
dimensions is governed by the Langevin equations for the 
rod center-of-mass position f = (3D ■ [Fu — W(r, </>) + f ] 
and for the rod orientation <j) = (3D r [M — d < pV{Y, <j>) + r] , 
respectively, where dots denote time derivatives and 
B^ 1 = ksT is the thermal energy. The rod's short time 
diffusion tensor D = -D||(u eg) u) + D±(I — u Cg> u) is 
given in terms of the short time longitudinal (-D||) and 
transverse (D±) translational diffusion constants, with 
u = (cos 4>, sin (/)), I the unit tensor and Cg> a dyadic prod- 
uct. D r is the short time rotational diffusion constant. 
Fu is a constant effective internal force that represents 
the propulsion mechanism responsible for the determin- 
istic motion in the rod orientation, and M is a constant 
effective internal or external torque yielding the deter- 
ministic circular motion (see the sketch in Fig. [I]). V(r, (j>) 
is an external confining potential, f and t are the zero 
mean Gaussian white noise random force and random 
torque originating from the solvent, respectively. Their 
variances are given by /|| (*') = 2S(t - t r )/(f3 2 D ll ), 
f±(t)f±(t') = 2S(t - t')/(/3 2 D ± ), and r(t)r(f) = 2S(t - 
t')/((3 2 D r ), where /||, f± are the components of f paral- 
lel and perpendicular to u, respectively. The bars over 
the quantities denote a noise average. We remark that 
for an active self-propelled particle, F and M are ef- 
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fective net forces that could be determined in the bulk 
from the forward and angular velocities F — \f]/(/3Dn) 

and M = \(f>\/(j3D r ), respectively, but are not necessar- 
ily directly connected to the internal propulsion mecha- 
nism [181 ]. 

At first we consider the free circle swimmer, i.e., we 
set V^r, (f>) = 0. In the limit of zero temperature, the rod 
center of mass would describe a perfect circle of radius 
R = (D\\F)/(D r M), with the circular frequency u> 
SI). \ I , For finite temperature all moments of r and 4> 
can be calculated exactly. The first and second moments 
of <p(t) are simply given by (f> = 4>$ + cot and A</> 2 = 

(tut) 2 + 2D r t, where 4>o — <H* = 0), and 
> run ad infinitum. The first two moments 
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where we let 
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with A = (5D\\Fj{D 2 + uj 2 ), u = (cos 0o , sin </) ) , 

= (— sin (/) , cos </>o )i u = (cos (f>, sin <f>), and u = 
(— sin0, cos <fi), i.e., Ar describes a spira mirabilis. 

We consider a very thin rod of length L, where 
D r /D\\ = 3/(2L 2 ), Dj_ = D||/2. We will denote all times 
in units of tb = L 2 / D\\ , lengths in units of L, and ener- 
gies in units of Different regimes are distinguished 
in terms of the dimcnsionlcss quantities D t /lo and (3FL. 
The latter determines whether the rod's erratic motion is 
dominated by the kicks of the solvent particles or by the 
self-propulsion. The former is the ratio of the ballistic 
over the random turning rate. In Fig. [T] we show Ar for 
different internal torques M and a typical trajectory of 
the rod position during two complete turns. In the second 
inset of Fig. Q] we display Ar 2 , which shows deterministic 
behavior for t < 1 / D r while for large times the swimmer 
moves in a random fashion according to Ar 2 oc t. 

Next, we introduce a confining, integrated segment- 
wall power-law potential in the x direction, V(x, <j)) = 
J^dlv[x'(l)} + kx with v(x') = (PL^HL/xT + 
[L/(L X — x')] n }, where L x is the channel width, n = 24 
is a large exponent, and x'(l) is the x position of the 
rod segment at contour length I (see the right inset of 
Fig. [2]). In case the solvent is confined as well, hydrody- 
namic interactions between the particle and the wall lead 
in principle to an x-dependent diffusion tensor [19( , which 
is ignored in our model. An additional gravitational force 
in the x direction [17| of strength k will be applied later, 
but we focus first on the symmetric case k — 0. At 
zero temperature, for a not too large ratio M/LF and 
under appropriate initial conditions (r o ,0o); the- tilted 
swimmer performs a steady-state sliding motion along 
either of the two walls with a constant x-position close 




FIG. 1: (Color online) Trajectories of the mean position r of 
the self-propelling rod for fixed f3FL = 10, f3M = 0.2, 1, 5, 25 
(ro = O,0o = 0). Left inset: the mean-square displacement 
Ar 2 for the same force and torques, but also for f3FL = 0, 
j3M = (lowermost curve). Right inset: a typical trajectory 
of the rod for f3FL = 25, /3M = 10, for times < t < r B . 
Lower right inset: Sketch of the self-propelled circle swimmer. 



to the wall and with a constant angle <j) determined by 
the steady-state conditions x = 0, and = 0, respec- 
tively. Without loss of generality, we consider the case 
M > 0, i.e., the rod rotates counterclockwise, such that 
it slides upwards along the left wall (see the sketch in 
Fig. [2]). In the limit of hard walls (n — * oo), the two so- 
lutions to the set of steady-state equations can be given 
explicitly as x s / u = L(l — 1/2 cos0 s / M ), (i.e., the front 
rod tip sits on the wall), and cos 2 <fi s / u = [1 — 2(M / LF) 2 ^f 
y/l - 8(M/LF) 2 }/[2 + 2{M/LF) 2 }, cos0 s/tl < 0, where 
the minus sign corresponds to the stable (0 S ) and the 
plus sign to the unstable {(j> u ) solution. Clearly, for 
2V2M/LF > 1 there is no solution to the steady-state 
conditions, but the rod keeps on rotating. For large ex- 
ponents n, the asymptotic steady-state velocity in the y 
direction is given by v y ~ D^F sin <f> s /(l + cos 2 <j> s ). 

The sliding mode is also present at finite tempera- 
ture. However, by thermal fluctuations the rod eventu- 
ally leaves the wall and reaches the opposite wall under an 
appropriate angle for the respective sliding mode in the 
opposite y direction, which we refer to as "flipping." Con- 
sequently, the circle swimmer moves diffusively according 
to Ar 2 ~ 2Di,t, with Dl the long-time transnational dif- 
fusion coefficient. This picture is clearly confirmed by 
Brownian dynamics computer simulations, averaged over 
1000 independent simulation runs, as shown in Fig. [21 

For large (3FL, large /3M, and a channel width of the 
order of the circle radius (L x < R), the average time the 
swimmer spends in its stable mode on either of the walls 
is large as compared to the duration of a flip. Thus, the 
swimmer effectively performs a one-dimensional random 
walk with a typical step length a ~ v y /j, where 7 is 
the flipping rate. This random walk leads to a long-time 
diffusion coefficient of Dl — v y/li which we display as 
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FIG. 2: (Color online) Mean-square displacement Ar 2 in con- 
finement {/3FL = 60, L x = 8L) without gravity and with zero 
torque (black), without gravity and with finite torque (red), 
and with torque and gravity (blue). The left inset displays 
the rod sliding along the walls. The right inset shows the 
confining potential without (left) and with (right) gravity. 




FIG. 3: (Color online) Long-time diffusion coefficient Dl as a 
function of the torque M for the swimmer in the bulk and in 
confinement for (3FL — 60 and L x = 6, 8, 10L. (a) Computer 
simulation, (b) rate theory. 



a function of internal torque M for different wall-wall 
separations L x in Fig. [3] 

It is clearly seen from the simulations [Fig. [21a)] that 
the diffusion in the channel is strongly enhanced as com- 
pared to the diffusivity of the free swimmer. In par- 
ticular, this strong enhancement is already observed for 
M = 0, as the narrow walls constantly align the rod in 
the y direction. However, the diffusion eventually slows 
with increasing wall- wall separation L x . For intermedi- 
ate M/LF « 0.15, diffusion is enhanced even further — in 
the simulations [Fig. [3Ja)] by an order of magnitude — 
displaying a much smaller dependence on L x . This non- 
monotonic behavior of Z?£ as a function of M is due to 
the stability of the sliding mode. 

To understand the nontrivial interplay of F, M, and 
L x in more detail, we identified from the simulations 
three different paths, (a), (b), and (c), dominating the 
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FIG. 4: (Color online) Top panel: the paths governing the 
flipping rate: (a) turning in the direction of the torque, (b) 
turning against the direction of the torque, (c) three-stage 
event — first, turning in the direction of the torque and then 
turning against it. Bottom panel: The rates of the three 
different paths as a function of M/LF for f3FL = 60, L x = 
8L. 



flipping rate 7. They all describe the transition from a 
stable mode at the left wall (4> s ,x s ) to another at the 
right wall (0 S + n, L x — x s ) due to fluctuations in the 
rod orientation (f>, whereas the translational motion just 
follows the internal force F and the confining potential 
V(r, (/)) [20j. These three different paths are sketched in 
Fig. |4] and are described as follows: The rod can slip out 
of its stable sliding mode by fluctuating in the direction 
of the torque [path (a)] or by fluctuating against it [path 
(b)]. In path (a), detachment from the (left) wall, which 
amounts to overcoming a barrier in the torque/angle from 
(j> s to 4> u , most likely also leads to finding the stable mode 
on the other (right) wall (for L x < B). Path (b), how- 
ever, is only successful if the rod orientation is subject 
to strong and fast fluctuations which enable it to make 
a turn of an angle (— ir + 4> u — 4> s ) before reaching the 
other wall. This explains why for intermediate torques 
and small L x , another important three-stage path (c) 
is dominating. This path is initiated by a small fluc- 
tuation of the orientation against the direction of the 
torque, from cj) s toward 7r/2 on the (left) wall. In a sec- 
ond stage, the swimmer approaches the other (right) wall 
at a small, constant turning velocity (j), reaching it after 
only a short time due to its strong internal force. By the 
other (right) wall it is reoriented in an upward direction 
before, in a third stage, turning quickly in the direction 
of the torque such that it reaches the original (left) wall 
at an angle <f> u . The flipping rate is now given by the 
path integral 7 oc J D(f>exp(— (3S[r, </>]/4), keeping initial 
and final configurations of <fi and x appropriately fixed. 
Here, the Onsager-Machlup action is given by 5[r, 4>] = 
f™dt' \d t ,^{t') - M + d^V {v{t')^{t'))\ 2 0, [U, with 
t' — (3D r t the normalized time. Note that our system 
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is non-Hamiltonian due to the internal driving force and 
the translation-rotation coupling. Hence, the least ac- 
tion path cannot be found as the minimum energy path 
in some energy landscape, as vastly studied in the lit- 
erature [23l . |24| |. In contrast, we now construct a non- 
Hamiltonian rate theory by assuming that — in the limit 
of large forces (3FL — the flipping rate 7 is dominated 
by either of the three paths [i = (a), (b), (c)], identi- 
fied in the simulation. The respective minimum actions 
are given by S^r,,^], with <fii(t') minimizing the action 
subject to the constraints [<^(0) = (j> s ,Xi(0) — x s ] and 
[4>i(oo) = 4> s ± 7r,Xi(oo) = L x — x a ], where the plus sign 
corresponds to paths (a) and (c), and the minus sign to 
path (b). Paths (a) and (b) [(c)] are further constrained 
by the condition not to reach (to reach) the initial wall 
between the initial and the final stage. 

In order to calculate the associated actions for the dif- 
ferent paths, we divide the trajectories into parts where 
the front rod tip sits on the original (left) wall and into 
parts where the rod moves at a constant turning velocity 
<\>i in between the walls. The former parts can then be 

expressed as the barrier heights 4 f£ m d(/)\M — d<j,V\ [HI, 
with <j) m — (f> u for path (a) and cf> m = tt/2 for paths 
(b), (c), whereas the latter are simply given by \dt>(j>i — 
■^| 2 *max> ''max being the normalized time it takes to swim 
from one wall to the other (£ max is chosen to minimize 
the action). The individual rates are roughly given by 
7, w exp[— j3Si/4]/TB, where the kinetic prefactors are 
crudely approximated by 1/tb, and plotted in Fig. 2] 
Summation over the individual rates yields the long-time 
diffusion coefficient Dl ~ w y7 _1 > with 7 ~ Y^ili plot- 
ted as a function of M for different L x in Fig. [3£b). The 
rate theory reproduces clearly the L x dependence and the 
nonmonotonicity of Dl as a function of M and attributes 
it to different rates of the paths (a) and (c). Moreover, 
the maximum in Dl is predicted to be weakly dependent 
on L x in agreement with the simulations. However, the 
actual values of the rate theory differ from the simula- 
tion data due to the crude approximation made for the 
kinetic prefactors. 

Finally, we study the effect of an additional gravi- 
tational field in the x direction (k > 0), breaking the 
symmetry of the channel potential (see the right inset of 



Fig. [2]) . On average, the swimmer is now situated more 
on the left than on the right channel wall, such that the 
sliding mode becomes ballistic (see Fig. [5]). 

In conclusion, we have studied the dynamic behavior 
of a self-propelled Brownian rod performing circular mo- 
tion. In the bulk, the analytical solution reveals long- 
time diffusive behavior. In channel confinement, an effi- 
cient stable sliding mode was identified that strongly en- 
hances the long-time diffusion along the channel as ob- 
tained by computer simulation and a non-Hamiltonian 
rate theory. If the channel is asymmetric, the sliding 
mode leads to ballistic long-time motion. 

The sliding motion of circle swimmers can be verified 
in experiments with different set-ups: First, catalytically 
driven nanorods [lCj, [llj and self-propelled magnetic col- 
loidal rods confined to a microchannel lrjj| will exhibit 
sliding [25f . Second, confined bacteria [aH Hi an d sper- 
matozoa£I|, Q move in two dimensions along circles. In 
fact, the typical radius of the observed circular motion 
is in the range of 10 — 1000/im for spermatozoae [1,J5[ 
and of the order of 50^m for Escherichia coli bacteria 0] . 
Therefore, the radii are typically larger but comparable 
with the particle sizes. When these particles are exposed 
to microchannels of similar widths as the observed radii, 
as realized for the bacteria the predicted huge accel- 
eration behavior should be observed, as has already been 
seen in 3D Third, vibrated polar granular rods [|| 
with an additional left-right asymmetry perform circle 
motions. When placed into a slit geometry, a sliding ef- 
fect may be observed here as well. 

Accelerating the dynamics in the channel by tuning 
the torque may be exploited as a mechanism to separate 
a certain species out of a crowded solution of different 
active particles. If a microfluidic channel is connected 
to a bulk mixture, the species moving quickest along the 
channel will arrive first at the channel end and can ef- 
ficiently be removed. This might be more efficient than 
traditional separation techniques such as capillary elec- 
trophoresis [261 ]. 
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